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ABSTRACT 



Aims. We propose a reflection model of the time delays detected during an exceptionally bright, single flare in MCG-6- 
30-15. We consider a scenario in which the delays of the hard X-rays with respect to the soft X-rays are caused by the 
presence of the delayed reflection component. 

Methods. We employ a model of the flare, which is accompanied by reprocessed emission. We consider two geome- 
tries/thermal states of the reprocessing medium: a partially ionized accretion disk surface and a distribution of mag- 
netically confined, cold blobs. 

Results. The reprocessing by cold blobs predicts positive time delays and a saturation in the time delay - energy rela- 
tion, which is likely present in the data. The model requires a strong reflection component and relies on the apparent 
pivoting of the combined primary and reflected spectrum. The reflection by the ionized disk surface does not reproduce 
the observed delays. We discuss the relation between the two reflection scenarios and argue that they are both present 
in MCG-6-30-15. 
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1. Introduction 

It has been known for a long time that Active Galactic 
Nuclei (AGN) show strong non-periodic variability in the 
X-ray band on timescales ranging from years down to hours 
and minutes (e.g. Lawrence et al. 1987; McHardy & Czerny 
1987; Mushotzky et al. 1993; Vaughan et al. 2003a; Kaastra 
et al. 2004). 

The predominant power-law character of the X-ray 
emission of most Seyfert galaxies and quasars can be suc- 
cessfully modeled assuming Compton scattering of soft pho- 
tons by a hot plasma (e.g. Shapiro, Lightman & Eardley 
1976; Sunyaev & Titarchuk 1980). In the accretion disk 
scenario of AGN the soft photons are naturally provided 
by the disk emission whilst the plasma can be identified 
with a hot corona of a temperature ~ 10 9 K. The detection 
of the high energy cut-off in a number of radio quiet AGN, 
with values between 60 and 300 keV (e.g. Zdziarski et al. 
1995; Risaliti 2002; Deluit & Courvoisier 2003; Soldi et al. 
2005; Beckmann et al. 2005), supports the thermal origin 
of the X-ray emission in these sources. 

The power-law shape of the X-ray spectra is a conse- 
quence of successive Compton up-scatterings. This should 
lead to systematic time delays between the hard and the 
soft X-ray variations as the hard X-ray photons must un- 
dergo more scattering events. A search for such delays was 
performed for MCG-6-30-15 and other Seyfert galaxies but 
the results are not unique. 



Reliable measurements of X-ray time delays in AGN 
spectra require sufficiently long observation times. For a 
95 ksec observation of MCG-6-30-15 with XMM-Newton 
the quality of the data did not allow to confirm any spec- 
tral delays (Ponti et al. 2004). They were only reported 
after a follow-up 300 ksec observation when Vaughan et al. 
(2003b) detected delays by 200 sec on timescales of 10 4 sec. 
Frequency- dependent time delays were found in a number of 
sources (NGC 7469, Papadakis, Nandra & Kazanas 2001; 
NGC 4051, McHardy et al. 2004; MCG-6-30-15, Vaughan et 
al. 2003b; NGC 3783, Markowitz 2005; Ark 564, Arevalo et 
al. 2006), for a frequency range of 10~ 6 to 10~ 3 Hz. These 
delays turn out to be roughly proportional to the Fourier 
period and they increase with the energy span considered 
(max. between 0.2 to 12 keV) up to a few thousand seconds. 

The dependence of the delay time on the Fourier fre- 
quency cannot be accommodated within the frame of a sim- 
ple Comptonization model that assumes electron clouds of 
a uniform density and temperature. In such models the ob- 
served hard-to-soft time delay is independent of the Fourier 
frequency (Miyamoto et al. 1988). Therefore, two types of 
more sophisticated models were proposed: 



an extended Comptonizing region with the appropriate 
density and temperature stratification (extended spher- 
ical region, Kazanas et al. 1997; advection-dominated 
inflow-outflow model, Blandford & Begelman 1999; 
extended coronal emission model, Arevalo & Uttley 
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2006) 



• models of a local spectral evolution due to coherent 
changes in a single compact emission region (single 
coronal flare evolution, Poutanen & Fabian f 999, Zycki 
2003; inflow of a single blob, Bottcher & Liang f999, 
Zycki 2004) 

The local models predict a pivoting of the hard X-ray 
spectrum, and the pivoting around the appropriate fre- 
quency can explain the apparent time delays, as argued 
by Kording & Falcke (2004). 

The detection of a single gigantic flare in the lightcurve 
of MCG-6-30-15 (Ponti et al. 2004) offers a specific insight 
into the delay problem assuming that in this exceptional 
case we see a localized flare event (instead of a random peak 
in the net flux of large flare distributions). This assumption 
is inspired by the rather symmetric shape of the lightcurve 
during the flare event and by its short time scale of only 
~ 2200 sec. Ponti et al. (2004) detected spectral time delays 
during the flare period and found that the hard X-rays were 
lagging the soft band. They interpret this as a result of 
multiple photon scattering within the flare itself. 

Using the same data, we discuss in this paper an al- 
ternative possibility to explain the time delays. While the 
models listed above focus primarily on the Compton up- 
scattering in the hot plasma, we investigate a possible role 
of the reprocessed component for the observed time delays. 
Unlike the primary, the reprocessed radiation arises in a 
relatively cooler medium and thus has a different spectral 
shape. The spatial distance between the primary source and 
the reprocessing site causes a delay between the two com- 
ponents and a time evolution of the combined (i.e. primary 
plus reflection) spectrum. We will show that this leads to an 
apparent effect of energy-dependent time delays. We inves- 
tigate two types of reprocessing: a strongly ionized reflector 
and a weakly ionized, cold reflector. 

After recalling the analysis of the data (Sect. 2) we de- 
scribe a simple model (Sect. 3) to find general trends of 
the spectral time delays between various energy bands. We 
discuss their dependence on details of the model setup for 
the two reprocessing scenarios (Sect. 4). We then show how 
the spectral delays of the observed flare can be reproduced 
by our model (Sect. 5) and we discuss and summarize our 
results (Sect. 6). 
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Fig. 1. Top: the light-curve over 0.2-10 keV during the pe- 
riod of the strong flare in MCG-6-3 0-15. Bottom: ene rgy- 
dependent time delays computed by iPonti et ail (|2004h to- 
gether with a linear fit representing Comptonization inside 
the flare. 

period is shown in Fig. [1] (top). It was decomposed into six 
lightcurves in separate energy bands, and X-ray time de- 
lays between them were calculated from cross-correlation 
functions (Ponti et al. 2004). The delays increase with dif- 
ference in energy (Fig. [TJ bottom) . The photons at 10 keV 
lag by a few hundred seconds with respect to those below 
0.57 keV. 



2. Observation and data reduction 

The data were collected during the observation performed 
on June 11 th 2000, with the EPIC MOS1, MOS2 and pn 
cameras of XMM-Newton (Striider et al. 2001) operated in 
timing, full-frame, and small window modes. The source 
was observed for 95 ksec and the spectral data were first 
published in Wilms et al. (2001) whilst Ponti et al. (2004) 
performed a timing analysis. 

The data were re-reduced and screened using the SAS 
software v6.5.0. We only used the pn data as there was 
no significant pile-up present. The single and double pixel 
events (i.e. pattern <4) were selected. Source and back- 
ground counts were extracted from two circular regions of 
45" radius. The background events were taken from a re- 
gion far from the source but within the same chip. 

During the observation MCG-6-30-15 showed a bright 
flare lasting for ~ 2000 sec. The lightcurve over the flare 



3. The model 

We assume that the flare is created in a compact region 
at some distance from the reprocessing matter. The result- 
ing primary component J p partly shines directly toward 
a distant observer and partly gives rise to the Compton 
reflection/reprocessed component 7 r . The response of I T 
to spectral variations of I p is delayed by the light travel 
time between the flare source, the reprocessing sites, and 
the observer. This time delay depends on the geometry 
in the vicinity of the flare and also on the viewing angle. 
An extended reprocessing region should lead to some level 
of smearing of the original signal, i.e. the duration of the 
brightening of the reflected component will be longer than 
the duration of the flare itself. 

We introduce a simple mathematical parameterization 
of J p and 7 r as a function of energy E and time t. For 
both components we assume a power-law dependence in 
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Table 1. Energy bands and fixed parameters of the model 



Energy band 


Range |keV] 


Fixed parameter Value 


ASi 


0.2 - 


0.57 


E m [ n 


0.2 keV 


AE 2 


0.6 


- 2.2 


^max 


12 keV 


AE 3 


2.2 


- 2.6 


Qp 


0.9, 1.3 


AE 4 


2.6 


- 4.8 


to 


2200 sec 


AE 5 


4.8 


- 6.8 


Ti 


500 sec 


AE 6 


6.8 - 


12.0 







We specify the energy bands, AEi, with i = 1, 6, and 
construct light-curves 

Li(t) = J I ohs {E,t)dE. (6) 

To investigate the time delay between the spectral 
response in different energy bands we compute cross- 
correlation functions 



energy with indexes a p and a r . For I t a broken power- 
law is also possible. It involves two branches with indexes 
a Tl and a r2 that connect at the break energy E\, T . This 
parameterization neglects particular reprocessing features, 
such as the emission lines, but represents rather well the 
overall shape of the continuum for a neutral or partially 
ionized reflector. 

We imagine the time modulation of I p and I T as consec- 
utive events. They are modeled by Lorentzians, £ p (t) and 
C r (t) , which take the half- width parameter Tf and the time 
to marking the maximum of the primary emission. The pri- 
mary leads the reprocessed spectrum by a delay 5 giving a 
measure of the average light travel time between the flare 
source and the reprocessing sites. As there are many possi- 
ble light paths for the reprocessing the width of C T can be 
broadened by a factor b. Note that this broadening is also 
influenced by the viewing angle of the distant observer. The 
equations for I p and I r read as follows: 



I p (E,t) = £ p (t)E- a * 
I T (E,t) = NC r (t) x 

where 



(1) 



br 



for^ < E hl , 
E- a *z forE > E bT , 



4>(*) = 

Cr(t) = 



T 2 
1 [ 



b 2 T 2 



(i - to - S) 2 + b 2 T\ 



2 ' 



(2) 
(3) 



A factor N is included in I r to normalize the reflected 
component against the primary. We define N with respect 
to the ratio, K, between I p and I T integrated over time and 
energy: 



Nb 



K = 



-Ebr 

/ dE'E'- 



J dE"E"- a -2 



(4) 



J dEE- a ? 



with b = [J C T (t')dt'] I [J C p (t)dt] . For a given value of K, 
the normalization N is a function of the energy limits, the 
broadening, and the spectral indexes. A distant observer 
detects the sum 



Iob s (E,t)=I p (E,t)+I r (E,t). 



(5) 



/ Ht)L- 3 {t-T)dt 



(7) 



I L 2 (t')dt'xJ f L 2 (t")dt" 



where the denominator is a normalization factor and T max 
must be large enough for Fq CF (t) to converge. The global 
maximum of F^, cf (t) marks the time delay between the 
energy bands AE\ and A£j. 



4. Results 

Our model contains various parameters. So me of them 
are fix ed by the flare observation described in iPonti et alj 
(|2004t ). The total energy range, and the energy channels are 
specified as in the analysis: E^in = 0.2 keV, E max = 12 keV 
and the energy bands are listed in Table [TJ We adopt the 
time of the flare maximum at to — 2200 sec and the half 
width Tf = 500 sec. The slope of the primary emission 
is not known very accurately due to the complex X-ray 
spectrum of MCG-6-30-15. It c orrelates with the hard X- 
ray flux, and ris es up to ~ 1.3 (|Vaughan &; Edelsonll2001l : 
IShih et aT]|2002f ). We thus choose a p = 1.3 as the upper 
limit of the slope. As the lower limit we set a p = 0.9, 
which is consistent with Bepp o-SAX data suggesting slopes 
- 0.94 dljb ian et all l2002al ) , and also wit h the average 
slope in a large sample of radio-quiet AGN dNandra et all 
Il997t iGeorge et al.ll2odot iPiconcelli et al.ll2005fh The fixed 
parameters of the model are summarized in Table [T] 

4.1. Delays for a highly ionized reflector 

Highly ionized reflection is likely to form when a strong flare 
is localized not far from the accretion disk. An example of 
the reflected spectrum for such a case is shown in Fig. O 
The sp ectrum is calculated usin g the codes Titan and 
NOAR (pumont et al.ll2000L [20031) . The disk vertical struc- 
ture is determined with an im proved version of the code 
given in iRozanska et"ai] (|2002t) . We assume M = 10 7 ' M 0) 
(iMcHardv et all 120051) and q/M = 0.998 dFabian et all 
2002at iMiniutti et al.H20Tjl iFabian fc Vaughanll2003f) . The 
flare is located at the disk radius r = 18i? g . 

The reflected component has a typical spectral shape 
for reprocessing by highly ionized material: it is relatively 
hard above ~ 3 keV and steep (soft) at lower energies: 



I T (E,t) = N£r(t) x 



£-1.57 for E < 3 A keV) 
^+0.15 iQlE > 3-1 ke y_ 



(8) 



The relative normalization of the primary and the re- 
flected component is taken from the modeling shown in 
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1000 10000 
Energy [eV] 

Fig. 2. An example of ionized reflection: the primary 
F P (E) oc E~ 09 (dot-dashed) and the local reprocessing 
(solid) of the ionized reflector scenario. The dashed line 
denotes the power-law fit F X (E) oc E~ a * adopted for the 
delay modeling. Below 3.1 keV the slope a Tl — 1.57 and 
above a r , = —0.15. 



Fig. O We obtain K ~ 0.15 if the flare radiates isotropi- 
cally. With b = 2 this leads to a normalization of N ~ 0.022 
(JV ~ 0.025) for a p — 0.9 (a p — 1.3). However, large values 
of the equivalent width of the iron line (Wilms et al. 2001, 
Reynolds et al. 2004) indicate that the reflection may be 
strongly enhanced with respect to isotropic emission (e.g. 
Fabian et al. 2002a, Ballantyne et al. 2003, Taylor, Uttley & 
McHardy 2003). Therefore, we keep TV as a free parameter 
of the model. 

We vary N between 0.01 and 0.9 and set the intrinsic 
delay S to 500 sec, 700 sec, 1000 s and 1700 sec. The spec- 
tral time delays at all energy bands AEi are obtained and 
normalized so that they are measured with respect to the 
one at AEi. The same normal ization was applied in the 
data analysis (|Ponti et al.ll2004T ). The results for a broad- 
ening of b = 2 and for the two values of a p are shown in 
Fig. El 

Obviously, for the very low normalization values of 
N = 0.022 or N = 0.025 no significant time delays ap- 
pear. Noticeable delays only exist for N being higher by at 
least one order of magnitude. But even at higher TV these 
delay curve s are v ery different from the ones obtained by 
iPonti et all ((2004). The model mostly produces negative 
delays, while the observed ones are always positive. 

The behavior of the delay curves as a function of N 
and S is partly analogous in both cases of a p . The overall 
normalization of the curves increases with S until a max- 
imum at ~ 1000 sec is reached, then, for even higher 5, 
the normalization decreases again. It also rises with N but 
saturates when N goes beyond unity. 

There are also some important differences between two 
cases of a p . For a p — 0.9 the delays are negative at all en- 
ergies and the shape of the curves is concave over the whole 
energy range. The minimum delay occurs around 3 keV. In 
the case of a p = 1.3 the shortest delays appear at a slightly 
lower value. The curves rise more quickly at intermediate 
energies and their curvature becomes convex around 5 keV. 
At high energies, long positive delays are obtained. 



We vary the broadening parameter b while fixing the 
time-integrated energy output of the reprocessed compo- 
nent. Thus, for each value of b we find the corresponding 
N from Eqn. Q. We consider seven values for b and iV 
with Nb = 1.8. The resulting delay curves for a p — 0.9 and 
S = 1000 sec are shown in Fig. |4j For a given value of Nb 
a whole range of curves is obtained. However, the charac- 
teristic concave curvature and the minimum around 3 keV 
are always preserved. Conducting the same investigation 
for other values of Nb delivers similar results. Therefore, 
varying b does not help to overcome the discrepancy be- 
tween the modeled time delays and those observed during 
the flare in MCG-6-30-15. 




Energy [keV] 



Fig. 4. Time delay curves derived for various values of b in 
the scenario of a highly ionized reflector. The normalization 
Nb of the total energy released is kept constant at Nb — 1.8. 
The bottom curve (red, solid) represents the value b = 1.5. 
From bottom to top, the value increases in steps of 0.5 for 
each curve and ends up at b = 5 for the highest curve (black 
dots /double-dashes) . 

4.2. Delays for a weakly ionized/neutral reflector 

In our second scenario we investigate the time delays ob- 
tained for a weakly ionized reflector. This is a likely situ- 
ation when the flare is located far from the accretion disk 
and/or the cold plasma is magnetically confined and denser 
than the plasma at the surface of an irradiated standard 
disk. 

As a specific example of the reprocessed component 
we use the model of irradiated magnetically confined, 
cold clumps. The radiative transfer calculation s for such a 
medium were performed bv lKuncic et al.l (|19970 . An impor- 
tant effect in their modeling is that successive reprocessing 
can significantly lower the spectral energy in the soft X-ray 
range leading to very hard spectral slopes of J r . This hard- 
ening depends on the number density of the clouds in the 
medium. Here we adopt the parameterization 

I I = NC I (t)E-° A , (9) 

which is inferred from the case of 50 individual clouds on 
the line of sight. The index a r = 0.1 is taken from fitting 
the corresponding curve in Fig. 2 of iKuncic et al.1 (|1997f ) 
between 0.2 keV and 12 keV with a power-law. We com- 
pute delay curves for 5 set to 500 sec, 700 sec, 1000 sec, 
and 1700 sec. The values of the normalization N are taken 
between 0.01 and 0.9 and the broadening parameter is set 
to b = 2. The results are shown in Fig. [5] 
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Fig. 3. Energy-dependent time delays for a highly ionized reflector with two slopes of the primary spectrum. Left: 
a p = 0.9, right: a p — 1.3. From top to bottom the panels represent intrinsic delays S = 500 sec, 700 sec, 1000 sec, 
and 1700 sec. The different curves denote the following values for iV: 0.01 (red, solid), 0.05 (orange, dots), 0.1 (yellow, 
short dashes), 0.3 (green, long dashes), 0.5 (blue, dashes/dots), 0.7 (pink, dashes/double-dots), and 0.9 (black, double- 
dashes/dots). Note that the vertical scale is different for both columns. 



As for the highly ionized reflector the very low values of 
N = 0.022 and N = 0.025 lead to no significant time delays. 
But for higher N the curves look very different from the 
previous scenar io. They resemble a lot more the observed 
delays found bv lPonti et alj (|2004f ). Note that strong reflec- 
tion was also required to model the energy-dependent frac- 
tional variability amplitude in MCG-6-30-15 (Goosmann 
et al. 2006), and it is consistent with the large value of 
the iron line equivalent width measured by several authors 
(e.g. Wilms et al. 2001, Reynolds et al. 2004). Therefore, 
the high value of N may be justified. 

The following trends can be derived from the model: 
the overall normalization of the time delay curves increases 
with J, goes through a maximum around S = 1000 sec and 
then decreases again. All curves have positive slopes and 
the spectra evolve from softer to harder X-ray energies. 
The shape of the curves changes with N. Higher values of 
lead to longer delays until N passes unity and the curves 
saturate at a level, which scales with 6. For the steeper 
primary spectrum with a p — 1.3 the delay curves show 
the same general features and tendencies as for a p = 0.9. 
However, for a p — 1.3 the overall normalization of the time 



delay curves is larger and the curves are generally steeper 
in the soft X-ray range. 

We investigate the influence of b in a series of calcula- 
tions with a p = 0.9, Nb = 1.8, and S = 1000 sec. The re- 
sulting time delay curves are shown in Fig. [5] For increasing 
values of b, the delays become smaller. This can partly be 
explained by the lower normalization N required to keep 
the same value of Nb = 1.8. It was shown above that the 
delay curves level down for lower values of N. Another point 
is that for larger b the variations of I T start earlier in time. 
Thus, the variations of I p are less ahead of the variations of 
I T . As for the highly ionized reflector, conducting the same 
investigation for other values of Nb delivers similar results. 

5. Modeling the measured time delays of 
MCG-6-30-15 

As shown in Fig. [3l the highly ionized reflector only predicts 
time delays that are non-monotonic with energy. This is due 
to the concave reflection spectrum of the scenario. On the 
other hand, a cold reflector with a hard single power-law 
spectral shape gives a monotonic increase. Certain param- 
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Fig. 5. Energy-dependent time delays for a weakly ionized/neutral reflector for two slopes of the primary spectrum. 
Left: a p = 0.9, right: a p = 1.3. From top to bottom the panels represent intrinsic delays S = 500 sec, 700 sec, 1000 sec, 
and 1700 sec. The different curves denote the following values for N: 0.01 (red, solid), 0.05 (orange, dots), 0.1 (yellow, 
short dashes), 0.3 (green, long dashes), 0.5 (blue, dashes/dots), 0.7 (pink, dashes/double-dots), and 0.9 (black, double- 
dashes/dots). Note the change in vertical scale for the last two panels of the right column. 
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Fig. 6. Time delay curves derived for various values of the 
broadening b in the weakly ionized/neutral reflector sce- 
nario. The normalization Nb of the total energy released 
is kept constant at Nb = 1.8. The upper curve (red, solid) 
represents the value b — 1.5. From top to bottom, the value 
increases in steps of 0.5 for each curve and ends up at b = 5 
for the lowest curve (black dots/double-dashes). 



eter sets of the latter scenario reproduce the observed time 
delays for the flare of MCG-6-30-15 quite well. In Fig. [7] 



we show two examples of satisfactory data representation. 
The model reproduces the curvature seemingly present in 
the data. This may indicate that it catches the apparent 
signal propagation better than the fit by a straight line 
used in Ponti et al. (2004). The primary slope of the two 
solutions shown was set to a p = 1.3 assuming that for the 
high flux state of the flare the spectrum should steepen 
(|Vaughan fc Edelsonll200ll : IShih et al.ll2002ft . 

The solution is not unique. Partly, this is due to the 
large error bars attached to the data points in Fig. [7] There 
are three free parameters involved: the normalization N, 
the broadening b, and the intrinsic delay S. For higher accu- 
racy of this modeling, the absolute normalization K of the 
reprocessed power with respect to the primary must be de- 
termined in an independent way. Then the value of Nb can 
be computed from Eqn.3](which links N and b). Obtaining 
K for the specific case of the observed flare in MCG-6-30-15 
would require a very detailed spectral analysis. A decom- 
position into primary radiation and multiple reprocessing 
components using a sophisticated radiative transfer model 
would be helpful. However, the flare period only lasted for 
~ 2000 sec and the obtained spectra are not of a sufficient 
accuracy to allow such an investigation. 

Nevertheless, we are able to give a few robust conclu- 
sions. The model indicates that an intrinsic time delay of 
1000 sec is required. Lower or higher values of 5 do not de- 
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is comparable to this value. Therefore, the size of the 
flare source has to satisfy 



Fig. 7. Satisfying representation of the measured time de- 
lays (red diamonds with error bars) for the flare/clumps 
scenario. The two curves rely on the following parameteri- 
zation: 5 — 1000 sec, b = 2.5, and N = 1.1 (solid, black) as 
well as S = 1000 sec, b = 2.5, and N = 0.7 (dashed, blue). 



liver any satisfying fits. The delay of 1000 sec translates to a 
light travel distance of ~ 60 i? s , assu ming a black hole mass 
M = 3.3 x 10 6 M Q (jMcHardv et al.ll2005h . The broadening 
factor 6 = 2.5 suggests a significant spread in the distances 
between the source and the reprocessor which would agree 
with the multiple scattering expected in the clumpy sce- 
nario. The lower limit of the normalization is at N = 0.7, 
hence the reprocessing has to be strong. In the cold reflec- 
tion model this implies a high covering factor of the flare 
source. 



6. Discussion and summary 

We considered two scenarios in which the time delays of 
the hard to the soft X-rays are caused by the delayed re- 
flection component. The dependence of the delay on energy 
is caused by a difference in the spectral shape between the 
primary and the reprocessed spectrum. In the first scenario 
the reflector is ionized, the reflection component is concave, 
the predicted time delays are non-monotonic with energy, 
and the observed time delays are not reproduced. In the 
second scenario the reflection spectrum is close to a single 
hard power-law and the predicted time delays increase with 
energy. In this case, the model reproduces qualitatively the 
observed time delay curve. 



6.1. Comptonization delay versus reprocessing delay 

Ponti et al. (2004) argued that the observed delays are 
caused by Comptonization inside the flare source, whilst 
in our approach we neglect the internal delay and model 
the effect of the delayed reflection. To justify this method 
the flare has to be compact enough. The lower limit for the 
flare size can be estimated from the compactness param- 



eter (see e.g. Cavaliere fe Morrison! 119801 : ISvenssorj [l994: 



iTorricelli-Ciamponi et al.ll2005[ ): I < 12. During the 95 ksec 
XMM-Newton observation, MCG-6-30-15 was exceptionally 
faint, and its 2-10 keV flux translates into a luminosity of 
L = 2.6 x 10 42 erg s -1 . The peak luminosity of the flare 



D fl > 6 x 10 1 



2.6 x 10 42 erg s" 1 



(10) 



so the light travel time across the flare itself can be as short 
as ~ 200 sec and this provides a lower limit for the inter- 
nal delay. Unfortunately, we cannot estimate the number of 
scatterings within the source because the observational lim- 
its for the extension of the X-ray emission in MCG-6-30-15 
are not conclusive (Beppo-SAX data indicate a cut-off en- 
ergy a t 100-470 keV, depending on the model; iFabian et ail 
I2002aft . 

In our successful weakly ionized/neutral reflector sce- 
nario the overall shape of the delay curves is rather ro- 
bust against changes of the parameters. We find that in 
most cases the curve slightly flattens towards larger en- 
ergy separations. This flattening is a specific prediction of 
the model and it is encouraging that it is indicated by the 
data. Comptonization models on the other hand predict a 
linear relation in the 0.2-12 keV band, i.e. much below the 
characteristic energy of the Comptonizing electrons. 

6.2. A flare located above the accretion disk 

The irradiation of a standard accretion disk by strong 
flare emission was studied in numerou s papers (see e.g. 
i Navakshin et all 120001: iNavakshinl 12000): [Ballantvne et alj 
I2OOU iRozanska et al.H2002t ICollin et alil2003[ ). If the flare 
is located near the inner part of the disk the disk surface 
is strongly ionized. Relativistic modifications of the flare 
spectrum and Doppler shifts then explain the presence of 
a broadened iron Ka line (see e.g. lDovciak et aDl2004al rbl) . 
For MCG-6-30-15 the relativistically distorted line varies 
less than the continuum, which can be accounted for by 
the light-bending model (jMiniutti fc Fabianl [2004] ) . This 
model is also considered for explanation of the soft-excess 
(Crummy et al. 2006, Petrucci et al. 2006, Ponti et al. 2006) 
although being confronted to other approaches (Gierlihski 
& Done 2004, Chevallier et al. 2006). 

We illustrate the simplified geometry of a flare scenario 
in Fig. [H] (left). The radiation is assumed to be beamed or 
bent toward the accretion disk and released within a cone of 
a certain half-opening angle 9q. A larger 6q leads to a wider 
range of light travel times between the flare source and the 
disk, and hence the broadening of I x increases. By varying b 
in our model we thus change the physical picture assuming 
a different angular release of the primary radiation. 

The geometry of the light-bending model assumes the 
flare source to be close to the disk axis and at varying 
heights above the disk surface. This lamppost geometry is 
justified by the idea of shocks during the relativistic jet 
formation (Henri & Pelletier 1991, Martocchia et al. 2000, 
2002). It is interesting to note that during the 95 ksec ob- 
servation the detected line possibly became narrower with 
time (<t = 0.39 ± 0.20 for the detection 3000 sec after the 
flare peak, and a — 0.19 ± 0.12 for the detection 4000 sec 
after the flare peak, Ponti et al. 2004), as if coming from 
more distant parts of the disk. 

The delay model does not constrain the radial posi- 
tion of the flare but it infers that the flare source is ele- 
vated rather high above the disk. Assuming that the flare 
is triggered by magnetic reconnection, the height of the 
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flare source above the disk surface should be equal to the 
disk's pressure scale height. This is implied by the equipar- 
tition between thermal and magnetic pressure as the flux 
tubes rise into the corona by buoyancy. An analytical fit 
of the disk's pressure scale he ight was given by Eqn. (9) 
in ICzernv fc Goosmann! (|2004h . With M = 3.3 x 10 6 M 
and assuming a high-state accretion rate of 0.3 times the 
Eddington rate (jUttlev et all [2002). we obtain pressure 
scale heights lower than 4.5 R g for radii R with 6 i? g < R < 
600 R g . From our delay model we infer an intrinsic time de- 
lay <5, which corresponds to an elevation of the flare source 
of 60 R g . Far from the axis, it is difficult to imagine a mech- 
anism transporting the flux tubes so high up. Therefore, a 
lamppost geometry where the source is driven upwards by 
a central outflow would seem more likely. However, as we 
discuss below there is another difficulty with the model in 
wh ich the primary sourc e o ccurs at such a large height. 

IBallantvne fc Fabianl (|2001h and iBallantvne et al.l 
(|2003l ) invoke two reflectors to explain the shape of the 
iron Ka line in MCG-6-30-15. They presume highly 
ionized reflection close to the black hole and a distant 
colder reflector at ~ 70R g . In principle, this cold reflection 
may also form if the flare is elevated high above the 
disk surface as proposed in the lamp post model. The 
reflection is weakly ionized or neutral if the ionization 
parameter £ = L/(Aird 2 n) < 10. Herein d is the distance 
of the flare from the disk surface with density n. With 
L = 2.6 x 10 42 erg s _1 and the local density in the disk 
atmosphere n — 10 14 cm -3 we obtain the condition for 
an almost neutral reflection as d > 1.5 x 10 14 cm, which 
corresponds to a light travel time > 1 hour. This is more 
than the measured delay time during the observation and 
makes the scenario of a highly-elevated flare unlikely. 

6.3. The weakly ionized reflector 

The question arises where a weakly ionized/neutral reflec- 
tor can be located if the flare is not allowed to be far away 
from the reprocessor. A possible answer is a model of cold, 
clumpy reflectors: the existence of magnetically confined, 
cold clouds in the accretion flow was suggested by Guilbert 
&: Rees (1988), and explored further by Celotti, Fabian 
& Rees (1992) and Kuncic, Celotti & Rees (1997). The 
interest in this model was recently revived by the detec- 
tion of narrow redshifted absorption lines which suggest an 
inflow over a narrow range of radii (NGC 3516, Nandra 
et al. 1999; Mrk 509, Dadina et al. 2005; E 1821+643, 
Yaqoob & Setlemitsos 2005; Q0056-363, Matt et al. 2005; 
PG 12111+143, Reeves et al. 2005; Mrk 335, Longinotti 
et al. 2006). The magnetically confined clumps were also 
recently discussed in the context of the very high state in 
galactic sources (see e.g. Yuan et al. 2006 and references 
therein). In general, partial covering remains an interest- 
ing alternative to the flat disk reflection for a number of 
Narrow Line Seyfert 1 galaxies (Mrk 335, Longinotti et al. 
2006; 1H 0707-495, Fabian et al. 2002b). 

An illustration of a clumpy reflector model is shown in 
Fig. [H] (right). We include some light trajectories between 
the flare source, the various reprocessing clouds, and the 
distant observer. The broadening effect on the reflection 
signal depends on the exact geometrical arrangement. In 
general, a range of light travel times around an average 
value is expected. Very low values of our model parame- 
ter b correspond to a geometrical arrangement, where the 



reprocessing clouds are distributed symmetrically with re- 
spect to the flare source and to the observer. By increasing 
b our setup drifts away from such a symmetry. 

During the flare the iron Ka line complex was not de- 
tected. A strong enhancement of the line emission was seen 
only ~ 3000 sec after the flare (Ponti et al. 2004). On the 
other hand, considerable spectral variability was seen dur- 
ing the flare. It is therefore possible that while the flare went 
off the line formed in the clumpy medium and was strongly 
broadened, as suggested by the parameter b we obtain from 
the delay modeling. Then, the line could not be identified 
in the data. The line signature became apparent only later 
when the flare radiation was subsequently reprocessed by 
more distant and larger parts of the disk. 

Our delay model requires a significantly enhanced re- 
flection which indicates a large covering factor of the flare 
source with reprocessing clouds. This is consistent with the 
(generally) large value of the equivalent width of the iron 
line in MCG-6-30-15, particularly during the low flux peri- 
ods (Fabian et al. 2002a). Malzac et al. (2006) pointed out 
that strong reprocessing can be realized by reflection in a 
clumpy medium, which presents an alternative to the classi- 
cal light-bending models, where the flare source is elevated 
above the disk. 

6.4. Conclusion for a global picture 

The observed flare is an exceptional event. Our model does 
not require the magnetically confined, cold clouds to ex- 
ist everywhere in the accretion flow. We rather suggest 
that MCG-6-30-15 exhibits both, ionized reflection from 
elements of an accretion disk and colder reflection from 
confined clouds. The particular flare event we discuss here 
possibly occurred in a relatively distant region dominated 
by magnetic confinement and weak ionization, which cor- 
rectly reproduces the observed time delays. Eventually, the 
magnetic confinement is correlated with the existence of 
the reconnecting flux tubes producing the flare. 

Evidences of occasional single flares occurring at dis- 
tances of the order of ~ 10 i? g — 100 R g from the black 
hole are given by the signature of their orbital motion (e.g. 
Iwasawa et al. 2004, Turner et al. 2006). Closer in, the 
continuous disk may recover again due to condensation, 
as suggested by Liu et al. (2006) and Meyer et al. (2006). 
Therefore, our model is not in contradiction with general 
models of various phenomena in this source based on the 
processes taking place in the innermost region. 

The qualitative results we obtain suggest to continue 
this line of research: to perform more sophisticated model- 
ing of time delays and to look for suitable single flares in 
X-ray lightcurves of this object and of other AGN. It would 
be particularly important to see whether the departure of 
the time delay - energy bin relation from a straight line, 
well reproduced in our model, is a characteristic property 
of the phenomenon. 
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